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ABSTRACT 

We consider the formation of dwarf carbon stars via accretion from a carbon AGB companion in light of the new 
107 object sample of Downes et al. (2004). This sample is now large enough to allow good mass determination 
via comparison of a composite spectrum to theoretical atmospheric models. Carbon dwarfs of spectral type M are 
indeed main sequence M dwarfs with enhanced metallicity and carbon abundance. We also calculate the predicted 
abundance of both M and of F/G carbon dwarfs, and show that the latter should be falsifiable in the near future. 



1. INTRODUCTION 

A carbon star is a star with carbon to oxygen ratio (C/O) 
greater than one and has a distinct atmospheric chemistry 
that yields a spectrum easily distinguishable from a typical 
(C/O < 1) star. Observationally, it has proven much eas- 
ier to detect large and intrinsically luminous asymptotic gi- 
ant branch (AGB) stars with C/0> 1 than it is to detect the 
much fainter carbon dwarfs at the low-mass end of the main 
sequence. As a result, it was generally assumed that carbon 
dwarfs were exotic objects, and that most carbon stars were 
on the asymptotic giant branch. Only recently has it become 
apparent that main sequence carbon dwarfs comprise not only 
a sizable population, but perhaps even the majority of all car- 
bon stars (Green 1992; De Kool & Green 1995; Green 2000; 
Downes et al. 2004). Since low-mass main sequence stars do 
not produce carbon in their cores and there is less carbon than 
oxygen in the interstellar medium, it is natural to ask how these 
carbon dwarfs form. 

The first main sequence carbon star detected, G77-61 
(Dahn et al. 1977), is a radial velocity binary with a cool 
white dwarf companion (Dearborn et al. 1986), and the car- 
bon dwarf stars PG 0824+289 (Heberetal. 1993) and CBS 
311 (Liebert et al. 1994) are spectroscopic binaries with hot 
white dwarf companions. These authors concluded that car- 
bon dwarfs may have formed via accretion from a carbon AGB 
companion. De Kool and Green (1994) showed that such a 
model indeed produces a carbon dwarf population more numer- 
ous than that of carbon AGB stars. 

In this paper, we consider the remaining issues with this 
model in light of new observations by the Sloan Digital Sky 
Survey. What remains to be done in order to establish that 
dwarf carbon stars are binary companions of former AGB stars? 
First, while carbon dwarfs have spectral type consistent with 
the late main sequence, it has not yet been possible to precisely 
determine their masses. If, for example, a solar mass carbon 
dwarf had spectral type M due to unusual spectral lines, then 
simulations would have assumed the star to be fully convec- 
tive (i.e., the entire star has C/0> 1) while for a solar mass 
star C/O need only be > 1 in a thin stellar envelope. Second, 
while a few carbon dwarfs with white dwarf companions are 
known, it has not been possible to find a companion for every 
known carbon dwarf. Combined with the small sample size of 
observed binary systems, this is not yet compelling evidence 
for the production of carbon stars by binary mass transfer; af- 
ter all, 42% of stellar systems with an M star as their primary 
are binary (Fischer & Marcy 1992). Certainly additional obser- 



vations showing that carbon dwarfs indeed have white dwarf 
companions would be compelling evidence for the binary mass 
transfer model. We can also search for other predictions that 
might provide the strong confirmation necessary to prove this 
model correct. 

In this paper, we solve both of these problems and present 
compelling evidence for dwarf carbon stars as post-mass trans- 
fer binaries. In § 2 we analyze spectra produced by the Sloan 
Digital Sky Survey of carbon dwarf candidates and attempt to 
quantitatively determine the properties of this population. In § 3 
we consider the dynamics of a binary system involving an AGB 
and an M star and whether the endstate of such a system might 
be a dwarf carbon star, as well as whether such carbon dwarfs 
would be produced in the observed frequency. In § 4, we con- 
sider the new predictions that such a model makes with an eye 
towards verifiability. Finally, in § 5, we consider the additional 
science made possible by the confirmation of this theory. 

2. ANALYSIS OF SDSS SPECTRA 

Recent observations from the Sloan Digital Sky Survey 
(SDSS) have been able to observe main sequence stars as dim as 
23rd magnitude (Stoughton et al. 2002; Abazajian et al. 2003; 
Abazajian et al. 2004) and therefore include a large sample of 
main sequence dwarfs. Because a carbon dwarf does not pro- 
duce TiO in large quantities, carbon dwarfs are uniquely iden- 
tifiable using the five filters applied by SDSS. In addition to 
the theoretically well understood carbon AGB stars compris- 
ing approximately 10% of the total AGB population, SDSS 
also finds that 0.3-0.4% of main sequence M stars show C/0> 
1 (Schroeder 2003). In addition, the survey has taken low- 
resolution (69 km/s pixel) spectra of dozens of carbon dwarf 
candidates, including one sample of 103 dwarf carbon stars 
(Margon 2002; Downes et al. 2004). We determine the param- 
eters of these stars by matching their spectra to simulated ones. 

Our model atmospheres for the synthetic spectra are based on 
using the ATLAS code (Kurucz 1995) and are built on the four 
standard assumptions: hydrostatic equilibrium, plane-parallel 
geometry, radiative equilibrium and local thermodynamic equi- 
librium. We calculate the synthetic spectrum for a range of at- 
mospheric parameters (T e ff, logg) in two main domains. These 
correspond to a G-dwarf (solar analog) and an M-dwarf. For 
each domain the carbon abundance is increased until the C2 and 
CN bands appear in the optical region and match in strength the 
observed ones. 

For both G- and M-dwarfs, the optical spectra attain the fea- 
tures of a carbon star spectrum when the carbon to oxygen ra- 
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tio exceeds 1.05. Both spectra were compared to the compos- 
ite carbon M dwarf spectrum produced by Schroeder (2003). 
Preliminary results show that the composite spectrum does not 
match the simulated carbon G star spectrum, but does match the 
spectrum of the main sequence M dwarf with enhanced C/O ra- 
tio. 

The spectra of the observed objects were also inspected for 
atomic absorption lines: particularly prominent are the Na I and 
K I resonance doublets, as well as Ca I and Mg I lines, and the 
Ca II triplet. Matching these, as well as the Balmer Ha and 
H/3 lines, with the synthetic spectra allows a rough temperature 
separation (G vs. M dwarf) and luminosity separation (from 
pressure broadening). A more difficult question is the overall 
metal content of these stars; the low resolution and S/N of the 
spectra prevent a detailed analysis; a comparison of our spectra 
with the analysis of two extremely metal poor carbon dwarfs by 
Cohen et al. (2004) indicates that our objects are not metal poor 
as a sample. 

We can conclude that the simulated carbon M dwarf spectra 
are a very good match for the composite spectra. Therefore, 
the population flagged by SDSS is indeed composed of M stars 
with C/0> 1.05 in the stellar photosphere. 

3. CARBON DWARFS AS POST MASS-TRANSFER BINARIES 

Analysis of SDSS spectra of dwarf carbon stars as done in 
§ 2 shows that indeed at least 0.3-0.4% of M stars on the main 
sequence are carbon stars and that those stars show signifi- 
cant enhancement in their carbon content. Our goal is to ex- 
plain how carbon dwarfs can comprise such a sizeable fraction 
of low mass main sequence stars. When carbon dwarfs were 
considered exotic objects, a likely candidate for their forma- 
tion was accretion from material originating in the photosphere 
of carbon-enriched asymptotic giant branch (AGB) stars. It is 
therefore natural to consider whether a more sizeable popula- 
tion might still be formed by mass transfer from AGBs. 

The sole known class of stars for which C/0>1 is common 
in their photospheres are asymptotic giant branch stars. AGBs 
include a carbon-oxygen core surrounded by a helium-burning 
layer, and the convective zone reaches into the layer of helium 
burning. However, during a thermal pulse, material from the 
carbon layer of the core will be dredged up and enter the con- 
vective zone, and as a result the photosphere will be carbon 
enriched. Since the material is composed of 22% carbon and 
only 2% oxygen by mass (Groenewegen & de Jong 1993), the 
photosphere can eventually attain C/0> 1 . As a result, between 
5 and 15 percent of AGB stars are carbon stars. 

In addition, an AGB loses mass at an incredible rate, in some 
cases as high as 10" 4 M Q /yr. Over the AGB lifetime of 10 5 yr, 
as much as 1 -2M Q may be ejected from the AGB. As a result, de 
Kool and Green (1994) have proposed that dwarf carbon stars 
might be the end product of a binary system involving an AGB 
carbon star and an M dwarf. Now that we have established that 
at least a factor of ten enhancement of carbon in the M star is 
required and that according to SDSS approximately 0.3-0.4% 
of M stars show such enhancement, we can consider whether 
these numbers are consistent with such an explanation. 

In order for this model to produce carbon M dwarfs of the 
type and in the quantity observed, we must show two things. 
We must show that at the conclusion of the AGB phase of the 
primary star, the M dwarf will be left with a C/O ratio that 
matches the observations in the previous section. We must also 
show that the rate of such binary systems yields the eventual 



0.3-0.4% fractional abundance of carbon stars among main se- 
quence M dwarfs. 

As the primary star becomes an AGB, it will overflow its 
Roche lobe and the binary system will enter a common en- 
velope phase (Iben & Livio 1993). During the common enve- 
lope phase, there is a net outflow from the AGB of as much as 
IO^Mq/vt with a typical velocity on the order of 20 km/s. We 
take the Bondi-Hoyle approximation that all gas passing within 
a cylindrical radius 



will be accreted by the secondary star. In the case of an M dwarf 
of mass 0.30M Q , this yields r acc =1.3 AU. For a binary system 
with separation 1 AU, r acc would subtend an angle of 53 de- 
grees, so the secondary will accrete material within a solid an- 
gle of 2.5 steradians. Thus, in the Bondi-Hoyle approximation, 
the mass Md of the dwarf changes as 

Md ~ -0.20Ma.GB- ( 2 ) 

As discussed in § 2, carbon M dwarfs show evidence of a 
range of metallicities, but as a whole are not metal poor. This 
is not surprising, given both that the elemental abundances in 
the photosphere of a typical M dwarf are solar (1) and that their 
low intrinsic luminosity restricts the SDSS sample to observa- 
tions of nearby dwarfs. So, we can take the initial abundances 
in the secondary star photosphere as (Aslpund et al. 2004) 
log C/log H = -3.61 ± 0.05, log O/log H = -3.34 ± 0.05. 

Since it is typical for both stars in a binary system to have 
formed simultaneously, if we take solar composition for the M 
dwarf, we should take the same initial composition for the AGB 
progenitor. Sivarani et al. (2004) observed the metal-poor car- 
bon AGB star CS 29497-030 to have [Fe/H] = -2.8, [C/Fe] = 
+2.38, and [O/Fe] =+1.67. Ohnaka et al. (2000) (see also 
(Lambert et al. 1986)) find the abundances in two carbon stars 
TX Psc and V Aql with [Fe/H] ~ -0.5 and -0.7 as [C/H] = 
0.16 and 0.13 and [O/H] = -0.14 and -0.36, respectively. If 
the surface gravity is lower than solar in TX Psc, Ohnaka et 
al. calculate an additional slight enhancement in C/O. All of 
the above abundances are given in dex relative to Anders and 
Grevesse(1989). 

As expected there is a metallicity dependence in [C/H] and 
[O/H], but for a star with [Fe/H] = 0, we can estimate from 
the above that [C/H] ~ 0.75 and [O/H] ~ 0.35. Since the ob- 
served carbon AGB stars may be at different stages and both the 
C/O ratio and the mass loss rate increase with time, this should 
be a lower bound on both the C/H and C/O ratios. These ra- 
tios in material accreted by the M dwarf companion are likely 
considerably greater than these bounds. 

If the AGB indeed loses a total of 1M Q over its lifetime, 
then using the initial abundances for a main sequence dwarf 
with mass 0.30M Q , the final mass of the secondary will be 
M k, 0.50M Q . Its abundances will then be 

[C///] d ~0.45, (3) 

[O/H] d ~0.17. (4) 

This results in C/O > 1 .02 as our lower bound. We have not 
applied the surface gravity correction (which, if analogous to 
TX Psc, would result in C/O > 1 .05), but the uncertainties in 
[C/H], [O/H], and [Fe/H] for carbon AGB stars are much 
larger than this correction. 

As shown in § 2, the C/O ratio needs to be enhanced until 
C/O > 1, or around a factor of five, to produce a carbon star. So 
we see that this model there is enough material accreted from 
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the AGB and with the proper abundance to produce a carbon 
star of the type imaged by SDSS. A star with C/O very close to 
1 will be an S star, showing prominent ZrO lines. Because the 
atmospheric chemistry is so sensitive to the C/O ratio around 
1, C/O = 1.02 is large enough to produce a carbon star rather 
than an S star. However, it should be noted that if further in- 
vestigation into the elemental abundances of carbon AGB stars 
reduce our lower bound, this might present a problem for the 
mass transfer model of carbon M dwarf formation. 

Now, consider the rate at which M stars are the secondary 
in a binary system in which the primary will undergo an AGB 
phase. In this case, we are interested in systems composed of 
a secondary of mass M s ~ 0.1 -0.4M Q and a primary of mass 
1.5 < M p /Mq. So, such a system has q = M s /M p < 0.4, and 
a typical q will be between 0.05 and 0.2. There are consid- 
erable difficulties in determining the frequency of binary sys- 
tems with different mass ratios, particularly for small values of 
q because the corrections for detection biases can be larger than 
the uncorrected signal(cf. Latham 2002). Indeed, surveys have 
found that there is a peak, with the frequency decreasing with 
smaller q (Duquennoy & Mayor 1991), that there are two peaks 
(Goldberg et al. 2003), or that the distribution is close to flat 
(Mazeh et al. 1992). It is also possible that some of the differ- 
ences between these samples are a result of examining different 
ranges for the mass of the primary. 

For each of these samples, the frequency of systems with 
larger q is much easier to measure and has a smaller required 
correction for biases. There is still considerable disagreement 
as to the shape of the distribution function f(q), but all of these 
distributions are consistent with 

^f(q = 0.8)dq < f(q = Q.\)dq < 3f(q = 0.8)dq. (5) 

Duquennoy and Mayor (1991) find that 57% of systems with 
F or G primaries are binary, while for comparison, Fischer and 
Marcy (1992) find that 42% of systems with M dwarf primaries 
are binaries. Taking f(q) flat to within a factor of three as above, 
this means that 

N(GM systems) N(FM systems) 



• 0.086 



(6) 



N(G stars) N(F stars) 

Additionally, from the Salpeter initial mass function 
(Salpeter 1955), the number density of stars goes as M~ 1,35 . 
Therefore, there should be approximately 4 times as many 
systems with M star primaries as G star primaries. Thus, the 
fraction of M stars in binaries with an F or G mass companion 
is 

N(GM systems) 



N(M stars) 



'0.021 



(7) 



Finally, the maximum radius of an AGB is around 1 AU. So, our 
model requires a binary with a separation of order 1-2 AU. The 
period distribution from Duquennoy and Mayor shows a near- 
Gaussian distribution in log(period) centered at log(r) = 4.8 
with g = 2.3. So, requiring a separation between 0.5 AU and 
2 AU as an estimate of our requirements, we find that the final 
fraction of carbon M stars should be 
Af(carbon M) 



N(M) 



'0.0021, 



(8) 



where again, the uncertainty is at least a factor of three. 

Finally, between 5% and 15% of AGB stars show C/0> 1 
(Wallerstein & Knapp 1998). Current models of AGB stars in- 
deed show that the carbon content in the stellar atmosphere 
increases throughout the lifetime of the AGB with additional 



thermal pulses, as does the mass loss rate. If this is indicative 
of all AGB stars showing carbon but for only 5-15% of their 
lifetime, then the above rate should be our final estimate of the 
frequency of carbon stars among the M dwarf population. If 
not, then our estimate must be reduced accordingly. The SDSS 
measurement of the fraction of M dwarfs with C/0> 1 is 0.3- 
0.4% (Schroeder 2003). Therefore, if indeed the majority of 
AGB stars have a carbon phase, our model is consistent with 
the observed frequency. 

4. ADDITIONAL PREDICTIONS 

In the previous section, it was shown that the predicted fre- 
quency of carbon stars for our model is consistent with obser- 
vations by SDSS. However, there are still large uncertainties in 
not only the binary mass ratio distribution function f(q), but 
also in the initial mass function. An improvement in either of 
these would reduce the uncertainty in our final frequency of car- 
bon M dwarfs to one much smaller than a factor of three, but 
at this point the predicted frequency is not by itself conclusive 
evidence for this model of their formation. So let us consider 
other ways in which this model might be tested. 

Perhaps the most striking prediction of a model in which 
carbon dwarfs are post-mass transfer binaries deals with the 
residual primary object. This model predicts that most carbon 
dwarfs will have a white dwarf binary companion with a sepa- 
ration of order 1 AU. If the AGB progenitor was an F or G star, 
the white dwarf will have a mass of order 0.5M Q , and therefore 

v orb 



15 km/s. 



sin ; 



(9) 



Depending on the angle of inclination, it should be possible 
to resolve such systems using Doppler velocities. Indeed, the 
radial velocity binary G77-61 (Dearborn et al. 1986) and the 
spectroscopic binaries PG 0824+289 (Heberetal. 1993) and 
CBS 311 (Liebertet al. 1994) are known to be carbon dwarf- 
white dwarf systems. To this point, the limitations on this 
search have been primarily sample size. Downes et al. (2004) 
find 107 carbon dwarfs with 15.6 < r < 20.8. It should be 
straightforward to observe these objects and search for a binary 
companion. A larger sample of carbon dwarf-white dwarf sys- 
tems would be compelling evidence for the binary mass transfer 
model. 

Another test comes from considering the end state of a sys- 
tem containing an AGB progenitor and a secondary G star small 
enough that it has not yet left the main sequence. Again using 
the Bondi-Hoyle approximation, the larger mass of the G star 
results in 

M dwarf w -0.39AW (10) 

However, a G star with solar metallicity is not fully convec- 
tive, and only the thin stellar envelope (M ~ 0.02M Q ) must be 
polluted. Therefore, the final atmospheric abundances of car- 
bon and oxygen will be very close to those of the carbon AGB 
donor. Therefore, this model predicts a considerable carbon F 
and G star population in addition to the carbon M star popula- 
tion. 

The fraction of carbon G stars is calculated using the tech- 
niques in § 3. In this case, the initial mass function does not 
enter because we are looking at G-G binaries. Further, the un- 
certainty in f(q) in this range is not a factor of three, but rather 
closer to 50% as (Duquennoy & Mayor 1991) 

N(GG systems) 



N(G primaries) 



:0.06±0.03. 



(11) 
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This results in 

0.006 ±0.003. (12) 

Because f(q) is known with much greater precision for systems 
with F/G secondaries, a good estimate of this population would 
at provide a better test for the binary mass transfer model than 
the fraction of M stars with C/0> 1. 

Unfortunately, while it has been easy to distinguish carbon 
M stars in the SDSS color space because of the lack of TiO 
lines, there is no similarly distinctive feature of a carbon F/G 
star. However, in the Downes et al. sample of carbon stars, 
a small fraction do show Balmer and Ca II lines and are be- 
lieved to be F or G carbon stars. Approximately 20 more of 
these objects were found by searching for metal-poor F sub- 
dwarfs (Schroeder 2003; Knapp 2004). Preliminary analysis 
(Knapp 2004) suggests that F/G carbon stars comprise 0.5% of 
the total main sequence population. This fractional abundance 
is consistent with the prediction in (12), and if confirmed in the 
final analysis would be very strong confirmation of the post- 
mass transfer model. 

5. SUMMARY 

We have used new observational evidence from the Sloan 
Digital Sky Survey to re-examine the model that dwarf carbon 
stars are formed via mass transfer from a carbon AGB com- 
panion. Using the SDSS composite spectrum of these stars, we 
have shown that carbon stars of M spectral type indeed are stars 
with C/0> 1 and that their masses are consistent with main se- 
quence M stars. Despite large uncertainties, we show that the 
predicted fraction of M stars with C/0> 1 is consistent with 
the observed fraction. This agreement is a requirement of our 
model but not by itself sufficient proof of its veracity. 

Two additional predictions should prove testable in the near 
future. There are already hints that many carbon dwarfs may 
be in binaries with white dwarfs. Now that there is a sample 
of over 100 carbon M stars, it should be possible to search for 
Doppler velocities we predict may be as large as 15 km/s. Such 
a search could likely be completed within two years. The other 



would be a search for F/G carbon stars. Large uncertainties in 
the mass ratio distribution of binaries and the initial mass func- 
tion apply much more strongly to M stars than F/G stars. As 
a result, agreement between the theoretical and observed abun- 
dance of F/G carbon stars would be strong evidence for our 
model. While F/G carbon stars do not lie in quite as distinctive 
a region of the color space produced by the Sloan Digital Sky 
Survey, there is already a project underway to search for F/G 
carbon stars, and this project should be able to provide a fairly 
stringent test for this model of carbon star formation upon com- 
pletion. Preliminary results suggest that indeed this project will 
provide the necessary confirmation of the binary mass transfer 
model. 

If this model for carbon dwarf formation is confirmed, we 
should turn our attention to its implications. One intriguing 
set of objects are the two very low metallicity carbon dwarfs 
discussed in Cohen et al. (2004). According to this model, a 
metal-poor carbon M dwarf is the endstate of a binary system 
with a metal-poor AGB and a metal-poor M dwarf. A com- 
parison of the abundances of the carbon M dwarf to those of 
an M dwarf of similar metallicity would allow a calculation of 
the elemental abundances present when its companion was a 
carbon AGB. Such a calculation has already been done with a 
metal-poor turnoff star to find abundances of its former AGB 
companion (Sivarani et al. 2004), so it should be feasible us- 
ing a carbon M dwarf as well. Therefore, one of the implica- 
tions of this model is that metal-poor carbon M dwarfs offer 
a probe into the elemental abundances of the earliest massive 
stars while they were on the asymptotic giant branch. While 
the main sequence abundances of these stars are of greater inter- 
est, particularly as a probe of primordial nucleosynthesis, these 
cannot be obtained without a considerable improvement in our 
understanding of the changes in elemental abundance between 
the main sequence and the asymptotic giant branch. 
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